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Within the framework of generalized factorization of higher-twist contributions to semi-inclusive 
cross section of deeply inelastic scattering off a large nucleus, multiple parton scattering leads to an 
effective medium-modified fragmentation function and the corresponding medium-modified DGLAP 
evolution equations. We extend the study to include gluon multiple scattering and induced quark- 
antiquark production via gluon fusion . We numerically solve these medium-modified DGLAP 
(mDGLAP) evolution equations and study the scale (Q 2 ), energy (E), length (L) and jet transport 
parameter (q) dependence of the modified fragmentation functions for a jet propagating in a uniform 
( — ^ ' medium with finite length (a "brick" problem). We also discuss the concept of parton energy 

, loss within such mDGLAP evolution equations and its connection to the modified fragmentation 

functions. With a realistic Wood-Saxon nuclear geometry, we calculate the modified fragmentation 
functions and compare to experimental data of D1S off large nuclei. The extracted jet transport 
parameter at the center of a large nucleus is found to be qo = 0.024 ± 0.008 GeV 2 /fm. 

PACS numbers: Valid PACS appear here 
Keywords: 

I. INTRODUCTION 

^ n 

Jet quenching or the suppression of large transverse momentum spectra [1[ can be used as an effective probe of the 
properties of dense medium created in high-energy heavy-ion collisions. Because of multiple scattering and induced 
gluon bremsstrahlung, an energetic parton propagating in dense medium will lose a significant amount of energy 
[243 an d therefore soften its final fragmentation functions. These modified fragmentation functions will lead to the 
' suppression of large transverse momentum single hadron spectra[l, Q, photon-hadron correlations (l0l - fl5| both away- 

side [TH, [l7| and same-side dihadron correlations [l8| in high-energy heavy-ion collisions. Such proposed phenomena 
have indeed been observed in experiments fl9l - |23j at the Relativistic Heavy-ion Collider (RHIC). Phenomenological 
studies of the observed jet quenching phenomena at RHIC indicate a scenario of strong interaction between energetic 
partons and the hot medium with an extremely high initial parton density [lfjl |2~H - [27| . The same phenomena are also 
predicted in deeply inelastic scattering (DIS) off large nuclei when the struck quark propagates through the target 
nuclei [28] though the extracted parton density in cold nuclei is much smaller than that in the hot matter produced 
in the central Au + Au collisions at RHIC. 

Large transverse momentum hadrons in high-energy nucleon-nucleon collisions are produced through hard parton 
scattering with large transverse momentum transfer and the subsequent fragmentation of energetic partons into final 
hadrons. Because of the large transverse momentum transfer involved the single inclusive, dihadron and gamma- 
hadron cross section at large transverse momentum can be calculated within the collinear factorized parton model 
of perturbative QCD (pQCD) [29rl3ll ]. Since the initially produced partons are extremely virtual they will have to 
go through a cascade of vacuum gluon bremsstrahlung before the final hadronization at a scale Qq ~ 1 GeV 2 below 
which pQCD is no longer valid. In the leading logarithmic approximation, the transverse momenta of the emitted 
gluons are ordered and therefore each successive gluon emission will contribute to one power of a logarithmic factor 
ln(Q 2 /Q 2 ) which will counter the small a s (Q 2 ) for large momentum scale Q 2 . The resummation of these vacuum 
gluon bremsstrahlung will lead to a scale dependence of the jet fragmentation functions which can be described by 
the Dokshitzcr-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [32[ evolution equations in pQCD. 

In the presence of nuclear and hot QCD medium, the initially produced energetic partons will have to go through 
multiple scattering and induced gluon bremsstrahlung as they propagate through the medium and before hadroniza- 
tion. The induced gluon bremsstrahlung effectively reduces the leading parton's energy and softens the final hadron 
spectra or parton fragmentation functions. Within the framework of generalized factorization of higher-twist contri- 
bution from multiple parton scattering, one can cast the modification of the semi-inclusive hadron cross section due 
to multiple parton scattering and induced gluon bremsstrahlung in terms of the effective medium modified parton 
fragmentation functions 0, 0| ■ I n the same spirit of the vacuum gluon bremsstrahlung, a medium modified DGLAP 
(mDGLAP) equation was also derived in Refs. @, which in effect resums successive induced gluon bremsstrahlung 
due to multiple scattering. The form of the mDGLAP evolution equation is the same as in the vacuum except that 
the effective splitting functions will contain both vacuum bremsstrahlung and the medium induced parts which should 
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FIG. 1: Gluon radiation from a single scattering (a) and double scattering (b-d) in DIS off a nucleus. 

contain information about the properties of the medium as probed by the energetic parton. 

Resummation of medium induced gluon bremsstrahlung in terms of a Monte Carlo simulation of parton shower 
has been carried out [34| [H| within the framework of opacity expansion for induced gluon bremsstrahlung [|| , which 
is equivalent to solving the mDGLAP evolution equations. The mDGLAP evolution equations as derived in the 
higher-twist approach jfjj, 0] have also been recently studied numerically (36j for modified fragmentation function in 
DIS off nuclear targets as this work is being completed. 

In this paper, we will extend the mDGLAP evolution equations as derived in Refs. [f| to include both induced 
gluon bremsstrahlung and induced quark-antiquark pair production from gluon-medium interaction within the same 
higher-twist expansion approach. We will then solve the coupled mDGLAP evolution equations for gluon and quark 
fragmentation functions using a numerical method based on Rungc-Kutta iteration. We will first consider parton 
propagation in a uniform medium with finite length and a constant gluon density (so-called "brick" problem). We 
will study in detail how the modified fragmentation function from mDGLAP evolution equations depend on the 
initial parton energy, the momentum scale, medium length and the jet transport parameter (or gluon density) . Since 
the concept of parton energy loss becomes ambiguous in the process of successive gluon bremsstrahlung and quark- 
antiquark pair production, we will consider the singlet quark distribution within a propagating jet with the initial 
condition of a J-function 5(1 — z) and a given flavor. We will study the momentum fraction carried by the singlet 
quark distribution function of the fixed flavor and compare to the quark energy loss as calculated in the traditional 
picture of induced gluon bremsstrahlung. With a realistic Wood-Saxon geometry of cold nuclear matter, we then 
calculate the nuclear modified quark fragmentation functions from the mDGLAP evolution equations in the DIS off 
nuclei and compare to experimental data to extract the jet quenching parameter in cold nuclear matter. 

II. MODIFIED FRAGMENTATION FUNCTIONS 

In DIS off a nucleus target, the struck quark via photon- nucleon interaction is likely to scatter with another nucleon 
along its propagation path inside the nucleus. Such secondary quark nucleon scattering can also be accompanied by 
induced gluon bremsstrahlung as illustrated in Figs. QJb)-(d). The induced gluon bremsstrahlung effectively reduces 
the energy of the leading quark jet and leads to a medium modified quark fragmentation function. We will work in 
collision frame where the momentum of virtual photon, the momentum of target nucleus per nucleon and the radiated 
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gluon are, respectively, 

q= l-Q 2 /2q-,q-,d T ], 
p=\p + , 0,0 T ], 

l = [%/2zq-,zq-,h], (1) 

and the Bjorken variable is x B = Q 2 /2q~P + ■ In this paper, we will use the convention for light-cone momentum (p) 
and coordinate (y) variables 

P = [p + ,P~,Pt], p + =p +p z ; P~ = (po -Pz)/2, (2) 
2/ = [y+jT,^], y + =t + z; y- = (t-z)/2. (3) 

Within the framework of a generalized factorization of higher-twist contribution to the semi-inclusive cross section 
from multiple parton scattering, medium modification to the quark fragmentation function was obtained from the 
double scattering in Figs. [ITb)-(d), the interference between single and triple parton scattering and the corresponding 
virtual corrections @, 0, H3 , 

AD h q (z h ,Q*) = [ Q §^ t ^[A^ qg (z,x B ,x L ,%)D h q {z h /z,Q 2 ) 

JO L T Z7r Jz h z 

+A-f q ^ gq (z 1 x B ,x L ,l^)D g ^ h {z h / z, Q 2 )], (4) 

where D q (z,Q 2 ) and D g ^h(z,Q 2 ) are the fragmentation functions in vacuum and the modification to the splitting 
functions, 



A"f q ^ qg (z,x B ,x L ,l%) - 



..2 



T A (x B ,x L ) + 5(l-z)A^T A (x B ,l 2 T ) 



(1 



NJ A (x B ) 

A%^ gq (z,x B ,x L ,l%) = Aj q ^ qg (l - z,x B ,x L ,l%) (5) 

are obtained from the induced gluon bremsstrahhmg spectra and therefore are related to the twist-four nuclear quark- 
gluon correlation distribution, 

T q A g (x B ,x L ) = J ^dyidy^e^+^^'il-e-^y^il-e-^^-y^) 

^\{AWMl + F+{y^)F + %y^)i )<l {y-)\A)6{-y^)e{ y --y^). (6) 

The matrix element, 

A^T A (xJ 2 T ) = f dz-^[2T A {x,x L )\ z=l -(l + z 2 )T A {x,x L )} , (7) 
Jo i- — z 

in the second term (which is proportional to a 5(1 - z) function) of the modification to the splitting function comes 
from virtual corrections to the multiple parton scattering cross section. This term can be constructed from the 
momentum sum rule (or momentum conservation) for the modified fragmentation function [(| 0, H3] , 

J dzzAD h q {z,Q 2 ) = 0. (8) 

The quark distribution function f^{x B ) of the nucleus is defined as 



= I ^e^ + ^i(A|^(0) 7 + ^( 2/ -)|A). 



With this modification to the fragmentation functions from single induced gluon emission, one can calculate the 
modified fragmentation functions by including the vacuum fragmentation functions which satisfy the vacuum DGLAP 
evolution equations, 

D h q {z, Q 2 ) = D h q (z, Q 2 ) + AD q (z, Q 2 ). (10) 
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If one defines the parton energy loss as the energy carried away by the radiated gluon, one obtains 

AE , a v a s f dll f , . , ,o . 

— = (Az g ) = — J -p- J dzzA-y q ^ gq (z,l T ) 



2 / a* t A C A l + (l-z) 2 T q A g (x B ,x L ) 



Assuming factorization of the two parton correlation distribution T qg in terms of quark and gluon distributions in 
nucleons (io . l4lj , one can express the twist-four nuclear quark-gluon correlation distribution, 



\ dy p N (y)[l - cos(x L p + y )] 



2%a s T~(xb,xl) 2ira s i A 
~N C f q A (x B ) = ~N~ C " 

x [{x L G N {x L ) + c(x L )[xG N {x)] x ~ ] 

dy~[l - cos(x L p + y~)]{q F (x L ,y) + c{x L )q F (0,y)\ 



(12) 



in terms of the generalized jet transport parameter in nuclear medium, 



qii(x L ,y) = ^ N 2 S _ * ' Pn{v)xlG n (x l ), (13) 

which in the limit of ij, = is also the transverse momentum broadening per unit distance for a parton in the 
^-representation of color. Here c(xl) = f q (xB + /^(xb)- The momentum fraction xl dependence in the 
generalized transport parameter reflects the energy transfer to the medium parton during the inelastic scattering with 
the propagating parton. Such recoil of the medium parton is in fact the cause for elastic energy loss 42j . We will limit 
the study to gluon bremsstrahlung and therefore will not consider the x l dependence of the generalized jet transport 
parameter q(xL,y) ~ q(0,y) = q(y). Under this approximation, 



2na s T^(x B ,x L ) 
N c L A (x B ) 



dy q F (y)4 sin 2 (x L p + y /2) . (14) 



Since Xl = oo for z = and both the quark J a [xb + x£) and gluon distributions xlGn(xl) vanish at xl = oo, 
the matrix elements related to the virtual correction in the modified fragmentation function becomes 

27ra s A^T A (x,e T ) f 1 ^ , 1 + z 



where 

,2* 1 



1 + Z f 

cb-j-— / dy-q F (y)4sm 2 (x LP + y- /2) , (15) 



<i T ) = 2 ( 1 - v 1 ~ 2l T/p +c r) ; 1 t < 2pV (i6) 

are the restrictions imposed on the limits of z and l T integrations by requiring xl < 1. Such limits are also imposed 
in the calculation of the modified fragmentation function. The radiated parton energy loss from single secondary 
scattering is then 

IT = Ca £ / dv ~ [ ^ f~* dz[l + {l-zY\ qF {y)Asm\x L p+y-/2). (17) 



III. MODIFIED DGLAP EVOLUTION EQUATIONS 



To take into account of multiple induced gluon emissions, one can follow the resummation of gluon bremsstrahlung 
in vacuum and assume that multiple medium induced bremsstrahlung can be resummed in the same way and one can 
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obtain the mDGLAP evolution equations for the modified fragmentation functions 0, 0] , 



dD h q {z h ,iJ?) a s (^) 



d In /i 2 



dD h g {z h ,y? 
d In fi 2 



2tt 



2tt 



1 dz 

, z 



1 dz 
, z 



lq^qg(z,p, ) D q{ — >M ) 



2n, 
9=1 



ff (*,/* a )^(^,/* a ) 



+ 7™(z,/. 2 )^(^,/i 2 ) 



(18) 



(19) 



where the modified splitting functions arc given by the sum of the vacuum ones 7 a _>.£ )c (z) and the medium modification 

A^ a ^ bc (z,l T ), 



la^bc{z, It) = la^bc(z) + Ma^bc(z, l T ) 



(20) 



with Aj a ^b c (z, l T ) for q — > qg given in Eq. ([5]). In the above mDGLAP evolution equations we deliberately used /x 2 
to denote the evolving scale whose maximum value in DIS is Q 2 

In addition to the first mDGALP evolution equation for the modified quark fragmentation function in Eq. (|18[) as 
derived in Refs. 0, 0] in DIS off nuclei, one has to consider multiple scattering and induced gluon bremsstrahlung 
for a gluon jet in order to complete the mDGLAP evolution equation for medium modified gluon fragmentation 
function in Eq. (TIT))) and the corresponding modified splitting functions. In the study of contributions from quark- 
quark (antiquark) double scattering to the modified quark fragmentation functions within the same framework of 
generalized factorization of multiple parton scattering [43|, one can in fact relate the effective splitting functions to 
the corresponding parton-parton scattering amplitudes. From gluon-gluon scattering matrix elements, one can obtain 
the following medium modification to the splitting functions (see Appendix 3 in Ref. [43| ) due to double scattering 
between a gluon jet and medium gluons, 



Aj g ^ qq ( Z , l T ) 



A lg^gg(zJ T ) 



w T [z2 + {l - z)2] 



l-^z(l-z) 



x / dy q F (y) 4 sin 2 (x L p + y /2) 
2C A N C (l-z+ z 2 ) 3 



z(l - z)+ 



IrpC F 



-5(1 -z)^ 



dy <?F(2/)4sm {x L p + y /2) 



A g (l T ) 



dy' 



l-e 



dz < n f [z 2 + (1 - z) 2 ] 



l-^z(l-s) 



2C A N C (1 -z + z 2 ) 3 



C f 



z(l-z) 



qF (y)4: sin (x LP y /2) , 



(21) 



(22) 



(23) 



where we have similarly assumed that the parton correlation distribution vanishes for xl > 1 which imposes limit on 
the integration over e(l T ) < z < 1 — e(l T ). The above calculations of the modification to the splitting functions are 
more complete since they include non-leading terms when 1 — z — > 0. For mDGLAP evolution equation for the quark 
fragmentation function, one can also compute the modification to the splittin g fu nction for q — ► qg from the complete 
matrix element of quark-gluon Compton scattering ((see Appendix 3 in Ref. [43|). 



A lq^qg(zJ T ) 



1 



C A 



z(l 



+ C F (l-z)(l + z 2 ) 



A q (l 



iK 1 t) 



dy 



(!-*)+ 

' q F (y)4: sin 2 (x L p + y-/2) - 5(1 - z) 
z(l + z 2 ) 



A q (l 2 ) 



dz 



C A 



+ C F (l-z){l + z z ) 



xq F (y)4: sin 2 (x L p + y /2) 



(24) 



(25) 
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Note that we have expressed the medium modifications to the splitting functions in mDGLAP evolution equations 
for both quark and gluon fragmentation functions in terms of the quark jet transport parameter qp and leave the 
color factors explicitly in the z-dcpendent parts. There is no single overall color factor for each effective splitting 
functions because they involve more than one channel of quark-gluon or gluon-gluon scattering which have different 
color factors. However, in the limit of soft gluon radiation, z — > 1 (here z is the momentum fraction carried by the 
leading parton), the modifications are dominated by leading term 1/(1 — z). In this case, 

A 7 „(z, I 2 ) « -^A 7g ^ qg (z, l 2 T ), for z->l, (26) 

and N C /C F = 2N 2 /(N 2 - 1) = 9/4. 

We can also similarly consider the quark-quark (antiquark) scattering and quark-antiquark annihilation processes 
for the secondary scattering. These processes will be responsible for flavor changing (conversion) in the parton 
propagation. They generally involve quark or antiquark density distributions of the medium and are also suppressed 
by a factor of l\jQ 2 [H| ■ We will neglect these processes for now and focus on the most dominant gluon radiation 
in the mDGLAP equations. 

In the following we will numerically solve the mDGLAP evolution equations for the modified fragmentation functions 
using a modified HOPPET (Higher Order Perturbative Parton Evolution Toolkit) 0] in LO. HOPPET is a Fortran95 
package with the GNU Public License that carries out the vacuum DGLAP evolution in z-space using Runge-Kutta 
method with a given initial condition D a (z,Q^). We replace the normal LO splitting functions in HOPPET with 
the modified splitting functions in each step of Runge-Kutta iteration to solve the mDGLAP evolution for modified 
fragmentation functions. 

Since all the medium modified splitting functions in mDGLAP evolution equations are proportional to the path 
integral of the jet transport parameter (?_f(j/) along the parton propagation length, the modified fragmentation func- 
tions will be determined completely by the spatial profile (and time evolution in the case of heavy-ion collisions) of 
the qF{y) which characterizes the properties of the medium that a jet probes. Therefore, we will study in the rest 
of this paper the dependence of the modified fragmentation functions on the profile of the jet transport parameter, 
its value, length of the medium, in addition to the energy (E) and scale (Q 2 ) dependence. By setting <?f(2/) = 0, the 
mDGLAP evolution equations will become those in vacuum and the corresponding vacuum fragmentation functions 
and their Q 2 dependence will be recovered. 

IV. THE "BRICK" PROBLEM 

We will consider first the simplest profile of the medium characterized by a constant value of q F = qo for a uniform 
medium with finite length L, assuming that a jet is produced at y = via DIS and propagates through the medium and 
finally hadronizes outside the medium. Such a "brick" problem is illustrative for understanding the main features of 
the modified fragmentation function from mDGLAP evolution equations. When solving mDGLAP evolution equations 
in this paper, we will use their LO form with rif = 3 number of quark flavors and a running strong coupling constant 
a s (Q 2 ) with Aq C d = 214 MeV. 

A. Initial conditions 

To solve the mDGLAP evolution equations one also has to specify the fragmentation functions at an initial scale 
Qq as the initial condition. Such initial conditions in principle should be different in medium and vacuum. Since 
perturbative QCD cannot be applied below the initial scale Qq the initial conditions are normally supplied by exper- 
imental measurements which are currently not available for medium modified fragmentation functions. Without such 
information, we will have to resort to a model of the medium modified fragmentation function at the initial scale Qq. 
Many studies [34], [3(| with the modified DGLAP evolution approach have simply assumed the initial condition at Qq 
in medium as the same in vacuum. This implies that partons below scale Qq will not interact with the medium and 
therefore suffer no energy loss. This is apparently counterintuitive in theory. Phenomcnologically, such an assumption 
also means that all medium modification comes from processes above scale Qq and therefore has a very strong Q 2 
dependence, which contradicts the experimental data as we will see later. To take into account medium modification 
to the fragmentation functions below the initial scale Qq, we will assume in this study, 

D a (Ql) = D a {Ql) + AD a (Q 2 ) , a = g,q,q. (27) 

where -D a (Qo) is the vacuum fragmentation function and AZ) a (Qo) will be generated purely from medium via the 
mDGLAP equations (|T5)) starting at fj, 2 = 0. In practice, we will simply set the effective splitting functions to contain 
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FIG. 2: (color online) Numerical violation of the momentum sum rule for the fragmentation function of a u quark as a function 
of the momentum scale Q 2 and for different values of the jet transport parameter qo- 



only the medium-induced parts, ^f a ^ c (z, lj) = A7 ^(, c (z, l 2 ^), freeze the running strong coupling constant at /i 2 = Qq 
and use the vacuum fragmentation functions as initial condition at /j 2 = in the calculation of AD^Qq). The value 
of Qq will be set at Qq = 1 GeV 2 in this study unless otherwise stated. 



B. Momentum sum rule 



To test the modified HOPPET method and estimate errors of the numerical solution to the mDGLAP evolution 
equations, we first check the following momentum sum rules 

[ dzzD*(z,Q 2 ) = l, (28) 
h Jo 

which should be satisfied for all parton fragmentation functions at any scale Q 2 . For this test, we assume the following 
initial fragmentation functions, 

D h q {z, Ql) = D h q (z, Ql) = D h g (z, Ql) = 5(1 - z). (29) 

which corresponds to a parton-hadron duality in vacuum at scale Q 2 with just one hadron species. Medium modifi- 
cation of the initial condition at Qq will be given according to the prescription in Sec. IIV Al 

Shown in Fig. [5J is the numerical violation of the above momentum sum rule for u quark fragmentation functions 

e= 1- f dzzD*(z,Q 2 ), (30) 
Jo 

in medium with different values of jet transport parameter q~Q. The sum rule for the fragmentation function in vacuum 
(<Zo = 0) is almost perfect and is also the case for small values of jet transport parameter qo < 0.1 GeV 2 /fm. The 
numerical error becomes larger as the value of 4q is increased in the "brick" medium. The HOPPET method uses 
a grid in z-space to numerically evaluate integration in z. This grid in z has a minimum value z m i n . Therefore the 
numerical integration over z is always truncated in z space below z m j„. Even though z m i n = e~ 30 in this study is 
a very small number but not exactly zero. So, the momentum fraction contained in the z < z m i n region causes the 
numerical violation of the momentum sum rule. Such violation should increase with Q 2 and qo because the mDGLAP 
evolution puts more hadrons (or partons) in the small z < z m i n region at larger Q 2 or for larger values of go- However, 
the sum rules in medium are still satisfied with an error better than 10 percent for most of the range of Q 2 and qo as 
shown in Fig. [2j 
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FIG. 3: (color online) Gluon distribution in the parton shower of a quark jet 
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FIG. 4: (color online) Sea quark distribution in the parton shower of a quark jet 



1.0 

z 



C. Shower parton distributions and energy loss 

Before we turn our attention to the medium modified fragmentation functions, it is useful to study the shower 
parton distributions generated by a propagating parton which can provide useful information about the underlying 
parton production and energy loss before the final hadronization outside the medium. To this end, we again assume 
parton-hadron duality for three species of "hadrons" , h = q,q 7 g and take initial conditions, 

D a a {z,Ql) = S(l-z);D b / a (z,Q 2 ) = (a,b = q,q,g), 

which simply mean that partons stop branching in vacuum at scale Qq. A parton with a large scale Q 2 > Qq, 
however, will undergo both vacuum and medium induced gluon brcmsstrahlung and pair production as described 
by the mDGLAP evolution equations in Eqs. (JTHJ) and (JTHJ). The original parton will then generate other species of 
partons within its parton shower. The initial condition for medium modification will be given by the prescription in 
Sec. I IV A[ which implies that jet-medium interaction continues to generate parton shower below scale Qq. 

Shown in Figs. [3] and @] are gluon D^(z,Q 2 ) and antiquark (or sea quark) distribution D^(z,Q 2 ) within a quark 
jet for different values of jet transport parameter qo, medium length L and momentum scale Q 2 . We have set Qq = 1 
GeV 2 and initial parton energy E = 100 GeV. Yields of gluons and sea quarks within a quark jet increase with 
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FIG. 5: (color online) Quark (antiquark) distribution in the parton shower of a gluon jet 



Q 2 for both vacuum (left panels) and medium modified (middle and right panels) fragmentation functions and the 
spectra also soften as more gluon brcmsstrahlung and pair production take place during the evolution with Q 2 . For 
small Q 2 the modified shower parton distributions become harder at large z due to medium-induced pair production 
of quark- ant i- quark and gluons (in the g — > gg splitting). However, with increasing Q 2 the modified fragmentation 
functions again soften and the yields increase with go and L due to the additional gluon bremsstrahlung and pair 
production induced by jet quenching during the propagation of the quark jet in medium. The sea quark distributions 
are generally softer than gluon distributions within a quark jet because quark-antiquark pairs are produced through 
gluon fission in vacuum and fusion in medium. 

Similar behavior is also seen in quark(antiquark) distribution D^(z, Q 2 ) within a gluon jet as shown in Fig. [5] For 
small Q 2 = 2 GeV 2 , one notes that the distribution is quite flat, reflecting the shape of the splitting function ; y g ^ q g(z) 
when the distribution is dominated by pair creation from the leading gluon whose initial distribution is peaked at 
3 = 1. As one increases qo, Q 2 and L, soft gluons are also generated which in turn produce soft quark-antiquark pairs. 
As a consequence, the final quark distribution within a gluon jet will become softer as the gluon distribution D^(z, Q 2 ) 
and sea quark distribution D^(z, Q 2 ) within a quark jet in Figs. OandSJ The gluon distribution within a gluon jet on 
another hand contains both the initial gluon and produced gluons from vacuum and medium induced bremsstrahlung. 
As one can see from Fig. [51 it indeed develops a peak at z ~ due to vacuum and induced brcmsstrahlung in addition 
to the peak at z = 1 of its initial distribution as one increases the values of %, L and Q 2 . Since one cannot separate 
initial and produced gluons, the concept of parton energy loss becomes ambiguous for a gluon jet in this picture of 
successive bremsstrahlung in medium. 

The case for a quark jet is different because net quark number is a conserved quantity. Even though the quark 
distribution contains both the initial leading or valence quark and the sea quarks from pair production, the sea quark 
distribution should be the same as the antiquark distribution and therefore can be subtracted. For the purpose of 
studying attenuation of the leading or valence parton due to vacuum and medium induced gluon bremsstrahlung we 
define the valence quark distribution as 



Shown in Fig. [7] is the valence quark distribution from the mDGLAP evolution equations. Because of gluon 
brcmsstrahlung and pair production during its evolution and propagation, the leading or valence quark distribu- 
tion in both vacuum (left panel) and medium (middle and right panel) gradually softens from its initial 8(1 — z) form 
as one increases the momentum scale Q 2 . Like other fragmentation functions, significant change of the valence quark 
distribution occurs already due to vacuum gluon brcmsstrahlung and pair production when it evolves from the initial 
Qq. As one increases jet quenching parameter qg and the medium length L, the valence quark distribution further 
softens. 

Such softening of the valence quark distribution can be characterized by the total fractional momentum change 



D v q (z,Q 2 ) = D«(z,Q 2 )-Dl(z,Q 2 ). 



(31) 
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FIG. 6: (color online) Gluon distribution in the parton shower of a gluon jet. 
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FIG. 7: (color online) Valence quark distribution in the parton shower of a quark jet 
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The net energy loss due to medium induced gluon bremsstrahhmg and pair production is then 

^ = AE - E AE " = fdzz [Dl(z,Qi)-Dt(z,Q>) . 



(32) 
(33) 

(34) 



Shown in Fig. ([5J is the medium induced energy loss as a function of the medium length L and initial jet energy E 
for different values of qo- Because of non-Abelian LPM interference, parton energy loss induced by multiple scattering 
as seen in the left panel clearly shows a quadratic dependence on the medium length 0]. However, for fixed initial 
parton energy E and jet transport parameter q , such quadratic length dependence gives away to a linear dependence 
for large values of medium length and the induced energy loss even saturates as AE/ E approaches to 1 because energy 
conservation imposed on each step of the mDGLAP evolution. 
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FIG. 8: (color online) Energy and length dependence of medium induced energy loss of a valence quark from mDGLAP 
evolution. 
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FIG. 9: (color online) Energy and length dependence of radiated energy loss of a quark from single gluon bremsstrahlung 
induced by multiple scattering. 



The behavior of such parton energy loss as defined by Eq. (|34|) with the modified fragmentation functions from 
mDGLAP evolution is quantitatively and sometime qualitatively different from the averaged radiative energy loss 
via a single gluon emission induced by multiple scattering. The averaged energy loss via single gluon emission is 
usually identified with the energy carried by the radiated gluon [Eq. (fl7|) ] as shown in Fig. |9] Even though energy and 
momentum is conserved in the gluon bremsstrahlung process as described by the higher-twist formalism, energy loss 
in this case of single gluon emission involves energy loss per emission multiplied by the average number of scatterings 
which could become a large number. Such a calculation apparently becomes invalid when AE/E > 1. This is an 
indication of the need for resummation in terms of mDGLAP evolution in which unitarity of multiple emission is 
ensured. 

In the calculation of energy loss in both the single gluon emission and the evolution of mDGLAP equations, we 
have imposed kinematic cut-off [Eq. (|16|) ] for the integration over fractional longitudinal and transverse momentum. 
This is the reason for the strong energy dependence of the energy loss for small values of E. The asymptotic energy 
dependence of the energy loss AE in both cases is weaker than linear though the results from mDGLAP evolution 
shows weaker E dependence than the single gluon emission. 
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FIG. 10: (color online) Modified it fragmentation function of a gluon (right), up, down (middle) and strange (left) quark jet 
that is produced at momentum scale Q 2 = 10 GeV 2 with initial energy E — 10 GeV and propagates through a uniform medium 
with length L = 5 fm and jet transport parameter go = (vacuum fragmentation function) (solid), 0.01 (dashed), 0.05 (dotted) 
and 0.1 GeV 2 /fm (dot-dashed). The initial condition to the mDGLAP evolution at Q% is given by HKN [S3] parameterization 
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FIG. 11: (color online) The modification factor RZ(z,Q 2 ) for the medium modified fragmentation function in Fig. 1101 



D. Modified Fragmentation Functions 



In the following, we will solve mDGLAP evolution equations with the initial conditions D^(z,Qq) given by the 
HKN parameterization of experimental measurements of fragmentation functions at Qq = 1 GeV 2 and obtain the 
corresponding medium modified fragmentation functions for hadron production. To quantify the medium modification 
of the fragmentation functions, we define the modification factor 

R ° M2) = §$Pr (35) 

as the ratio between modified and vacuum fragmentation functions. 

Shown in Figs. UOland ll H are the modified 7r + fragmentation functions for gluon (left panel), u, d (middle panel) and 
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FIG. 12: (color online) Modified it fragmentation function of a gluon (right), up, down (middle) and strange (left) quark jet 
with initial energy E = 10 GeV that propagates through a uniform medium with length L = 5 fm and jet transport parameter 
qo = 0.05 GeV 2 /fm. The momentum scale of jet production is Q 2 = 1.0 (solid), 2.5 (dashed) ,10 GeV 2 (dotted) and 100 GeV 
(dot-dashed). 

s (left panel) quark jets and the corresponding modification factors, respectively, for different values of jet transport 
parameter go in a uniform medium with length L = 5 fm. The initial jet energy is set to E = 10 GeV with virtuality 
Q 2 = 10 GeV 2 . The modified fragmentation functions are seen to be suppressed at medium and large momentum 
fraction z with increasing values of qo because of increased induced gluon bremsstrahlung and pair production. The 
produced soft gluons and sea quarks will also hadronize and contribute to the hadron spectra, leading to enhancement 
of medium modified fragmentation functions at small z. Since the medium modified gluon splitting function is about 
9/4 of a quark [Eq. (|26|) ] at large z, it radiates more soft gluons and consequently its modified fragmentation is more 
suppressed at intermedium and large z as compared to the quark fragmentation functions. 

At large z, the modification factor for gluon jets is seen to have a small bump for small values of q. This feature 
results from the competition between gluon bremsstrahlung and pair production in the evolution of a gluon jet. Both 
of two processes contribute to the vacuum evolution. In the medium, however, gluon bremsstrahlung suppresses the 
fragmentation function while pair production is relatively enhanced in particular at large z. If one switchs off the 
medium induced g — > gg and g — » qq channels, the bump disappears. This feature is similar to the shower parton 
distributions of a gluon jet as discussed in the last subsection. However, as one increases the value of q, these features 
at large z disappear as induced gluon bremsstrahlung further suppresses the gluon fragmentation functions. 

Since the modification to the effective splitting functions due to gluon bremsstrahlung induced by multiple scattering 
are power suppressed at large Q 2 , most of the the modification to the fragmentation functions should come from 
mDGLAP evolution at small Q 2 . Shown in Figs. [12] and [13] are the modified fragmentation functions and the 
corresponding modification factors for different Q 2 but for fixed initial energy E and value of qo =0.05 GeV 2 /fm. 
The fragmentation functions at Qq = 1 GeV 2 correspond to the medium modified initial conditions that we generate 
according to our assumed prescription in this paper. One can see in Fig. 1131 that the initial fragmentation functions 
are already suppressed at large z due to induced gluon emission and therefore most of the modification comes from 
mDGLAP evolution near or below the initial value of Q 2 = Qq. Medium modification above the initial scale is 
relatively small. This will give a weak Q 2 dependence of the medium modification of the fragmentation functions as 
we will show in the discussion of jet quenching in DIS. The modification, however, varies with the initial energy E 
significantly, as shown in Figs. [T4l and IT51 for a fixed value of qo = 0.05 GeV 2 /fm. Such variations reflect the energy 
dependence of the energy loss as shown in Fig. [5] Eventually, when the initial energy becomes infinitely large, the 
modification will become increasingly smaller. 
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FIG. 13: (color online) The modification factor RZ(z,Q 2 ) for the medium modified fragmentation function in Fig. 1121 




FIG. 14: (color online) Modified n fragmentation function of a gluon (right), up, down (middle) and strange (left) quark jet 
with initial energy E — 10 (solid), 20 (dashed) and 50 GeV (dotted) that propagates through a uniform medium with length 
L = 5 fm and jet transport parameter qo = 0.05 GeV 2 /fm. The momentum scale of jet production is Q 2 = E 2 . 



V. MODIFIED FRAGMENTATION FUNCTIONS IN DIS 



To calculate the nuclear modified fragmentation functions in semi-inclusive DIS off a nucleus we have to employ a 
more realistic form of the nuclear density distribution. We consider the initial quark jet produced at yo that travels 
along a direction with impact parameter b (see Fig. [If)] for illustration). We assume that the jet transport parameter 
along the quark jet trajectory is proportional to the nuclear density, 

~, ,n » PA(y,b) 

q{y ' b) = qo ^o-y (36) 

where pA(y,b) is the nuclear density distribution normalized as J dyd bpA(y,b) = A and go is defined to be the jet 
transport parameter at the center of the nucleus. We will use the three-parameter Wood-Saxon form of nuclear density 
distribution here. 
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FIG. 15: (color online) The modification factor RZ(z,Q 2 ) for the medium modified fragmentation function in Fig. 1141 




FIG. 16: Illustration of the path integration in DIS. See details in text. 



If we neglect the nuclear and impact parameter dependence of the nuclear quark distribution function, the photon- 
nucleon cross section that produces a quark at (yo,b) is proportional to the nuclear density distribution pA{yo,b). 
Then the averaged mFF should be 



D{z) = (D{z,y )) = - db 

A JO J -co 



dy Q D(z,y )p A (y ,b). 



(37) 



In order to calculate the D(yo, b) for a quark produced at location (yo, b), the path integral in the modified splitting 
functions should be replaced by the following 



dyg{y)[- 



PA (0,0) J y0 



dypA(y,b)[ 



(38) 



Shown in Fig. [T7] are the calculated nuclear modification factor of w , K and pip) fragmentation functions with 
different values of jet transport parameter qo — 0.016 — 0.032 GeV 2 as compared to the HERMES experimental data 
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[43 |. The experimental results are presented in terms of multiplicity ratio R^, which represents the ratio of the 
number of hadrons of type h produced per DIS event for a nuclear target of mass A to that for a deuterium target 
0: 

( Eejq f (x)D h f (z) \ / ( Ee 2 q f (x)D h f (z) 
{ Xe 2 q f (x) lA ) I \ Ve)q f {x) lD 

In the experimental data, each z bin has an averaged value of initial jet energy (E) = (u) and virtuality (Q 2 ) which 
we use in the calculation. Since the averaged values of Q 2 = 2.25 — 2.65 GeV 2 in HERMES data are quite small, the 
suppression factors are apparently sensitive to the value of the initial evolution scale Qq. In the calculation, we have 
used Qq = 1.0 GcV 2 . We used the HKN [5(J parameterizations as initial condition which extends to down to scale 
Q 2 = 1 GeV 2 . The agreement between the calculation from mDGLAP evolution and the HERMES experimental data 
[49j are quite well for three nuclear targets in the intermediate region of z values. At large values of z the agreement is 
not so well maybe due to other effects such as hadronic interaction [5l|, [52j that we have not taken into account. The 
calculated modification factors for protons in a large nucleus are also much smaller than the experimental data. This 
might be related to the non-perturbative baryon transport in hadronic processes [53] ] . We have so far neglected quark- 
anti-quark annihilation contribution to the mDGLAP evolution equations. These processes will affect the medium 
modification of anti-quarks and will likely improve the modification factor for anti-proton distribution. 

Shown in Fig. (|18j) are the calculated nuclear modification factors at fixed z as a function of the initial jet energy E 
as compared to the HERMES experimental data |4{|. Again, each bin of E has an average value of (z) and (Q 2 ) which 
we also use in the calculation. Similarly as illustrated in the uniform medium (a "brick"), the medium modification 
of the fragmentation function gradually disappears as the initial jet energy E increases. The agreement between our 
theory calculations and experimental data is generally good for a range of values of qo = 0.016 — 0.032 GeV 2 , except 
at lower energy where hadronic absorption might become important. From the combined fit we find the jet transport 
parameter at the center of large nuclei is qo w 0.024 ± 0.008 GcV 2 /fm. This value is consistent with the transverse 
momentum broadening of the Drell-Yan dilepton production in p + A collisions (54- 56] (Ag|,) = 0.016A 1 / 3 GeV 2 
which gives an averaged jet transport parameter qo = 0.018 GeV 2 /fm. 

As described in Sec. IIV A[ we have modeled the initial condition for fragmentation functions in the medium at 
Qo = 1 GeV 2 as different from the fragmentation functions in vacuum due to energy loss of parton nearly on shell. 
Therefore, most of the medium modification to the fragmentation functions come from mDGLAP evolution at low Q 2 
while contribution from high Q 2 region is power-suppressed. This will lead to a very weak Q 2 dependence as shown 
in Fig. 1191 where we compared the calculated modification factors for the fragmentation function with experimental 
data at fixed z and initial jet energy E but as a function of Q 2 . The calculated suppression factors are almost 
independent of Q 2 , consistent with the experimental data (49| . If one has chosen the initial condition for medium 
modified fragmentation functions at Qq as the same as the vacuum one, one would obtain a modification factor that 
has a too strong Q 2 dependence to be consistent with the experimental data. 

VI. SUMMARY 

In this paper, we have extended the modified DGLAP evolution equations to include induced gluon radiation for 
gluon jet and quark-anti-quark pair creation from gluon fusion (43j within the framework of generalized factorization 
for higher-twist contribution to multiple parton scattering. The effective parton splitting functions are proportional 
to a path integration of the jet transport parameter q over the propagation length. We then numerically solve the 
coupled mDGLAP equations for medium modified fragmentation functions for different static profile of medium and 
different values of the jet transport parameter q. For a "brick" profile, we have studied the shower parton distributions 
within a jet during its propagation through the medium and relate the change of momentum fraction carried by the 
valence quark as the effective energy loss of a propagating quark. We have systematically studied the energy, medium 
length and virtuality dependence of the medium modified fragmentation functions. Because of the inclusion of induced 
pair creation and different splitting functions for gluon and quark jets in the mDGLAP evolution, we found that the 
modification factors are different for gluon, valence quark and sea quark fragmentation functions. 

We also applied the mDGLAP evolution to quark propagation in the deeply inelastic scattering (DIS) of a large 
nucleus and found the calculated nuclear modification of the effective fragmentation functions in good agreement with 
experimental data in the intermediate z region. In modeling the initial condition for modified fragmentation functions, 
we have chosen to include medium induced radiation and parton energy loss below the initial scale Qq. In this case, 
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most of the medium modification comes from mDGLAP evolution in the low Q 2 region while large Q 2 contribution 
is power-suppressed. This leads to a weak Q 2 dependence of the medium modification of the fragmentation functions 
which is consistent with the experimental data in DIS. 
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FIG. 17: (color online) The modified multiplicity ratios as a function of z with different values of the jet transport parameter qo 
compared with the HERMES [4tJ data for Ne,Kr and Xe targets. For clear presentation the modification factors for different 
targets have been shifted vertically by some value(Kr by -0.05 and Xe by -0.2). 
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FIG. 18: (color online) The energy dependency of the nuclear modification factors with different values of the jet transport 
parameter qo compared with the HERMES [i^] data for Ne, Kr and Xe targets. For clear presentation the modification factors 
for different targets have been shifted vertically by some value(-RTr by -0.05 and Xe by -0.2). 
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FIG. 19: (color online) Comparison of the modified multiplicity ratios as a function of Q 2 at fixed value of z and jet energy E 
with the HERMES [431 data for Ne, Kr and Xe targets. For clear presentation the modification factors for different targets 
have been shifted vertically by some value(_Kr by -0.05 and Xe by -0.2). 



